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ABSTRACT 

The energy lost by a rotation-powered pulsar is carried by a relativistic flow 
containing a mixture of electromagnetic fields and particles. In the inner regions, 
this is thought to be a magnetically dominated, cold, electron-positron wind that 
is well described by the MHD equations. However, beyond a critical radius r^, 
the same particle, energy and momentum fluxes can be transported by a strong, 
transverse electromagnetic wave with superluminal phase speed. We analyze the 
nonlinear dispersion relation of these waves for linear and circular polarization, 
and find the dependence of on the mass-loading, magnetization and luminosity 
of the fiow, as well as on the net magnetic fiux. We show that, for most isolated 
pulsars, the wind lies well outside Tcr, and speculate that superluminal modes play 
an important role in the dissipation of electromagnetic energy into nonthermal 
particles at the termination shock. 

Subject headings: plasmas - waves - acceleration of particles - pulsars:general - stars: 
winds, outflows 
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Introduction 



Pulsar winds are believed to be launched as relativistic outflows consisting mainly of 
electrons and positrons. These winds are Poynting-fiux dominated when they are launched 
from the pulsar. Outside the light- cylinder (located at p = r/rLc = rVL/c= 1, with Vl the 
angular velocity of the pulsar) they move radially, carrying a frozen-in, nearly transverse 
magnetic field that is modulated by the pulsar rotation. They can be viewed as nonlinear 
waves, whose phase speed equals the radial flow speed (and is, therefore, sub luminal). The 



19901), in which the 



simplest analytic description of this flow is the "striped wind" ( ICoronitil . 
magnetic field is purely toroidal, has a magnitude proportional to 1/r, and reverses sign 
across a corrugated current sheet that separates the two magnetic hem ispheres. Numerical 
solutions of the force-free, oblique rotator problem (ISpitkovskyI 
support 



20061 ) ■ lend qualitative 



(Petri & Kirk 


2005; 


Petri 


2011) 



201 ll ) as well as the orbital modulation of the high anc 



very 



201 ih . 



high energy emission from gamma-ray binaries that contain pulsars (iPetri fc Dubus 



In the absence of dissipation, these subluminal waves propagate at constant speed 
and remain Poynting dominated, provided t heir density is high enough to apply the 



MHD approximation ( iKirk fc Mochol 



2011al ). However, basic dynamical considerations 



indicat e that the winds should not be Poyn ting dominated after they cross the termination 



shock ( jRees fc Gunn 



1974 



Begelman 



19981 ). and spectral and morphological modeling 



of pulsar wind nebulae suggest an even stricter upper limit on the Poynting dominance. 
How Poynting flux is converted to kinet ic energy flux is still not clear (for reviews see 



Gaensler fc Slane 



2006 



KirketaL 



20091 1. Reconnection in the current sheet, as proposed 
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BP 53 F-38041 Grenoble, France 
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bv ICoronitil (119901) and 



Kirk fc Skjaeraasen 



2003 



Lyubarsky 



Vlichell (1994), pro ceeds rather slowly ( iLyubarsky fc KirkI 12001 



20101) • In the case of an isolated pulsar such as the 
Crab, it fails to convert a significant fraction of the electromagnetic flux before the flow 
reaches the termination shock , unless the r ate of pair injection by the pulsar is much higher 



than conventionally assumed (lArons 



2011I ). 



Motivated by this finding. 



LyubarskyI (120031 ) suggested that the fields are dissipated 



in the termination shock itself. This is not permitted at an MHD shock in a strongly 
magnetized flow, where the compression ratio remains close to unity. Nevertheless, it 
is possible to imagine a viable scenario in which reconnection in the current sheet of a 
striped wind is driyen or triggered by the relatively weak compression of the MHD shock 



(^Petri fc Lyubarsky 



2007; 



exhibit such an effect (jS 



Lvubarskv fc Liverts 



than those suggested by 



AronsI (12011 



roni fc Spitkovskv 



20081 ). Recent 2D and 3D PIC simulations 



20111), albeit at even higher pair injection rates 



An alternative picture is offered by superluminal waves. Because the density in the 
wind falls with increasing radius, these modes, which are absent in the MHD description 



1975 



Melatos fc Melrose 



19961). 



can propagate outside of a certain critical radius ( lUsovl 
If this radius lies inside the termination shock, the possibility arises that the flow converts 
from a subluminal to a superluminal wave either spontaneously at some point upstream of 
the shock, or as an integral part of the shock itself. The latter possibility is an attractive 
scenario, since these modes are known t o damp much mo re rapidly than the subluminal 



1978 : 



Asseoetal 



1980l ). guaranteeing a relatively thin 



striped wind mode (iLee &: Lerche 
transition region. 

At low amplitudes, tra nsverse, superlu minal waves in an electron-positron plasma 



have simple properties (e.g.. 



Iwamoto 



19931 1. However, under pulsar conditions, the waves 



are strongly nonlinear. Their properties have been investigated in many papers, (for 
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example 


Kaw & Dawson 


1970: 


1984 




Melatos & Melrose 


1996 



Max & Perkins 


1971 


; 


Glemmow 


1974 




1977 




Asseo et al. 


Skjaeraasen et al. 


200^ 


)). Here we take a new look at 



the problem, and extend the known results for both circularly polarized waves, and 
linearly polarized waves with a non-zero (phase-averaged) magnetic field, using the 
two-fiuid (electron-positron) model. In particular, we formulate and solve the equations 
that determine into which modes a striped wind of given magnetization can convert. 
This amounts to solving the electrodynamic equivalent of the jump conditions across a 
magnetohydrodynamic shock front. As a result, we find the cut-off radius of these modes as 
a function of the striped wind parameters. We also show that the frequently used "strong 
wave limit" does not provide an adequate description of these modes. 



The paper is set out as follows: in Section [2] we review the basic equ ations for 



waves. Using the "homogeneous" or "H-frame" formulation introduced by 



Clemmowl (jl974j . 



ane 



19771 ). we compute the fiuxes of particles, energy and momentum carried as functions of 
the wave parameters. In Section [3l we use these relations to solve the jump conditions 
that must be fulfilled when the wind converts from a subluminal to a superluminal mode 
across a transition layer that is narrow compared to its distance from the pulsar. Section H] 
contains a discussion of these results and their relation to the problem of dissipation at the 



termination shock. 



2. Superluminal waves in an electron-positron plasma 



An analytical trea t ment of non-linear waves in a cold plasma was first undertaken 



by lAkhiezer fc PolovinI (1l956l ) who considered a one-component electron plasma in a 
background of fixed ions. In their approach, all wave quantities depend only on the phase 
(j) of the wave. For waves propagating in the x-direction, cf) = uj{t ~ x/v^), where is the 



phase velocity and u the frequency of the wave. Both subluminal and superluminal waves 
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(f^ < c and v^> c respectively) are possible, and the longitudinal and transverse modes are, 
in general, coupled. Pulsars are most likely surrounded by an electron-positron pair plasma, 
possibly loaded with a small number of protons, so that a single-fluid approximation is 
not valid. Also, the wave amplitudes near the pulsar are expected to be high enough to 
drive even the proton component of a plasma to relativistic velocity, suggesting that the 
rest mass of the fluid particles is unimportant. Consequently, rather than considering the 
general case of an electron-ion plasma, attention has been concentrated on the much simpler 
two-fluid electron-positron plasma. In such a plasma, both linearly and circularly polarized 
waves with pur ely transverse fields can propagate witho ut introducing fluctuations in the 



1975 



Kennel fc Pellat 



19761 ). It is these waves that we consider 



charge- density (lAsseo et al. 
in the following. 

Far from the pulsar (at p ^ 1), the magnetic field is twisted into a tightly- wound 
spiral form, and it is an excellent approximation to neglect the phase-averaged radial field 
component compared to the toroidal one. Since pulsar winds are radial, this means that the 
phase-averaged component of the magnetic field is transverse. Furthermore, in this region, 
the radius of curvature of the wavefront is much larger than one wavelength, so that the 
waves can be approximated as locally plane. 

In the following, we identify the (radial) propagation direction of the wave as the x-axis 
and the toroidal direction as the z-axis, so that = Q and = ^- The equations to be 
solved are the continuity equations and the equations of motion for each of the two particle 
species, along with Maxwell's equations. These are given in Appendix [XJ where it is shown 
that, from Maxwell's equations and making use of the assumption of purely transverse 
waves, the two species have the same density: 7+n+ = 7_n_ (a plus (minus) subscript 
denotes positrons (electrons)). Furthermore, the Lorentz factors are equal: 7+ = 7_ and 
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the other components of the dimensionless four-velocity of the two species are related by: 



To keep the notation simple, we henceforth drop the plus and minus subscripts and express 
everything in terms of the positron fluid quantities. 



For superluminal waves, 



Clemmowl (11974 , 



19771 ) showed the problem is significantly 



simplified if one performs the calculations in a frame of reference that propagates with 
a speed c^/v^ with respect to the lab. frame. In this frame, which we call the H-frame, 
the wave quantities do not depend on space but only on time, and the phase of the wave 
is just = Lot. An immediate consequence of Faraday's equation is that in this frame 
the magnetic field, in addition to being homogeneous, is also constant. In the following, 
unprimed symbols denote quantities in the H-frame; we use primed symbols for quantities 
in the lab. frame. 



2.1. Circular polarization 

The circularly polarized wave is described in the H-frame by three (phase-independent) 
parameters: the (constant) proper number density of each species n = hq, the component 
of the four- velocity in the direction of propagation Ux = u^o, and the magnitude of the 
four-velocity component perpendicular to the propagation direction ^/Uy + u^, which, for 
this wave mode, equals the dimensionless strength parameter denoted by a and defined in 
terms of the electric field amplitude Eq as 

a = e \Eo\ I {mcu) . (1) 

Consequently, the Lorentz factor of the fluids is also phase-independent: 7 = 7o = 



yTT"u|^~+~a^. The frequency u coincides with the (proper) plasma frequency, Up-. 



u = Up = \ . 2 
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In terms of these quantities, the particle flux density J and the (0,0), (0,a;) and {x,x) 

components of the stress energy tensor are given by (see Appendix |Bl equations IA37tiA39l 
and K40lE42D : 

J = no-foci—^] (3) 
V 7o / 

T°° = no^omc' (2^0 + (4) 

= no7omc^ {2um) (5) 



T'' = no7omc'[^ + -] (6) 
7o 70/ 



and are also phase- independent. 



2.2. Linear polarization 

The linearly polarized mode requires an additional parameter to fix the amplitude of 
the phase- averaged magnetic field. In the H-frame, the four parameters are: the proper 
number density tiq and the component u^o of the four velocity in the direction of propagation 
both taken at phase zero, (the point in the wave where the magnit ude of the electric field is 
maximal, and = Uy = 0), the nonlinearity parameter q, used by lKennel fc PellatI (Il976l ). 
and the ratio A of the (phase- independent) magnetic field to the electric field Eq at phase 
zero. Physically q is double the ratio at phase zero of the energy density in the fluids to 
that in the electric field, as measured in the H-frame: 

327rno7omc2 

1 = ^2 • (7) 
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Defining the normalized field variable y = E/Eq, one finds: 

Ux = M^o + 470(1 -Z/)A/g (8) 

7 = 70 + 270 (1 - y') /g (9) 

Uy = (10) 

= \/7^ - - 1 

= 2jo^/N{y)/q (11) 



where A^(?/) is a fourth-order polynomial, as follows from (|8]) and ([9]), and 70 = \/l + m^q- 
Periodic wave solutions exist provided N{y) has four distinct real roots yi..A, in which 
case y oscillates between the values 2/2 and y^ (assuming the ordering yi < y2 < ys < y^j- 
Because of the normalization of y, one has j/3 = 1 and 1^2! ^ 1- The strength parameter of 
the wave, as defined in ([T]) is 



dyi/v^W) (12) 



y2 



and the wave frequency is given by 
The phase- aver aged fluxes are: 



u = — —Up . [16) 



(J) = w(^) (14) 
/T-oo\ 2/0 , 47o(|/^ + A^) \ 

^T""") = no7omc ^27 H ) (15) 

(r°i> = no7omc^(2.x. + ^) (16) 

2/2w' 47o(y^+A2)\ 

T ) = ?2o7omc ( ^ ) (17) 

\ 7 Q / 



and the phase-averaged electric field is 



no7omc 
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Ampere's law enables these averages to be expressed as integrals over the normalized 
electric field y: 

dy Aiy)-i / JN{y) 

(Aiy)) = (19) 

(see Appendix [XI). 

If the phase-averaged magnetic field vanishes (A = 0) one finds y2 = —1, and the 
integrals in f[T^ - f[TH]) can be expressed in closed form in terms of elliptic integrals. For 
A 7^ 0, closed forms can be found to lowest order in an expansion in the small parameter q 



see 



Kennel fc Pellat 



19761 ). However, these are not adequate to describe the pulsar case, 
as we show in the following section. In the general case, one must resort to numerical 
integration, noting that the integrands have integrable singularities at each of the limits. 



2.3. Conserved quantities 

In order to identify those parts of the wind of a given pulsar in which superluminal 
waves can propagate, one has to associate the wave properties with wind quantities inferred 
from the observational data. A pulsar wind transports particles, energy and magnetic fiux, 
and exerts a ram-pressure on its surroundings. Across a stationary shock front the quantities 
J', T'°^ and T'^^ are conserved. Also, from Faraday's law, the transverse electric field 
component E' is conserved. Across a thin transition front where one wave-mode converts 
into another, and which is stationary on the oscillation timescale, the phase-averaged values 
of these quantities are conserved. 

However, it is the relative magnitudes of these quantities, rather than their absolute 
values, that determine the physics of wave propagation. We therefore introduce three 
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dimensionless quantities: 



V 



mc (J') 

mc{J') 



(20) 
(21) 
(22) 



mc{J') 

which are conserved across a transition layer and a re independen t of the density parameter 



no- /i is the mass-loading parameter introduced by 



Michel! ( 119691 ) ■ which corresponds to the 



Lorentz factor each particle would need if they were to carry the entire energy flux. 

For the superluminal wave modes discussed above, these quantities can be computed 
by Lorentz transforming their values in the H-frame into the pulsar, or lab. frame. A boost 
of (dimensionless) speed /3> = c/v^ in the negative x direction is required, for which explicit 
expressions are given in Appendix [Bl (We use the notation /3> and r> for the speed and 
associated Lorentz factor of the H-frame of the superluminal wave, and /3< and r< for the 
phase (or bulk) speed and associated Lorentz factor of the subluminal wave, both as seen 
from the lab. (pulsar) frame.) Applying these boosts, one obtains, for linearly polarized 
modes, expressions for /i, v and r], as functions of the four input parameters Uxo, q, A and /3> 
that must be evaluated by numerical integration. The case of circular polarization, where 
ry = 0, is much simpler, since fj, and z/ are given in closed form as functions of u^o, a and /3>. 

For the cold subluminal wave modes, such as the striped wind (in the absence of 
dissipation), fi and u are independent of r adius in the inner par ts of the wind where a 



single-fluid MHD description is adequate (iKirk &: Mochol 



estimated from standard pulsar models (e.g.. 



2011a 



Lyubarsky fc Kirk 



Their values can be 



200ll ). However, these 



two quantities are almost equal in highly relativistic winds, and it proves more convenient 
to use, instead of u, the magnetization parameter a, defined as the ratio of magnetic to 
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particle energy flux according to: 



V 



r<(i + a) 

r|(l + a)-(l + a/2) 



(23) 



(24) 



In the striped wind, the third parameter, r/, depends on the phase-averaged magnetic field 
value 



Without loss of generality, we may choose {B') > 0, so that < 6 < 1. Both h and r] are 
functions of latitude in the striped pulsar wind, since the spacing in phase of the current 
sheets varies. At latitudes greater than the inclination angle between the magnetic and 
rotation axes, conventionally denoted by a, the sheets vanish (zero spacing), whereas on 
the equator, they are equally spaced (separated in phase by tt). 

Thus, for given /i, a and 6, the jump conditions can be solved for three of the four 
input parameters Ma;o, 5', & and /3>, provided the fourth is specified. The wave frequency, 
normalized to cjp, can then be constructed from f[T3|) . 




(25) 



via the expression: 



a/3<r 



(26) 



3. 



Results 



3.1. Circular polarization 



For circular polarization, the jump conditions can be solved analytically ( iKirkI 120101 ). 
An example is plotted in Fig. [H which shows the refractive index ck' /u' = /3>, the strength 
parameter a and the component of the particle speed in the propagation direction, u^o, as 
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Fig. 1. — The circularly polarized modes that satisfy the jump conditions for /i = 12, a = 3 as 
a function of lab. frame frequency normalized to the proper plasma frequency. No solutions 
exist below a cut-off frequency that is close to w'/wp = 1. Two solutions, one with Uxo ~ fi 
and one with small u^o (denoted by solid and dashed lines, respectively) exist above the 
cut-off, up to oj'/up ~ 1.36, where a —t- for the Iow-m^^q mode. The modes have the same 
refractive index for given u'/up, which equals the speed /?> of the H-frame seen from the 
lab. frame. 
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a function of lab. frame frequency u' normalized to the (proper) plasma frequency Wp. In 
this example, we have chosen moderate values of the energy and momentum fluxes in the 
lab. frame: /i = 12, a = 3 (and, therefore, r< = 3), in order to display the structure of 
the solutions near the frequency cut-off. Since the wavenumber k vanishes in the H-frame, 
the refractive index is trivially obtained from the Lorentz transformation of the wave 
number and frequency: u' = T^co and k' = f3^T~^u/c, (this also applies, of course, to 
the case of linear polarization) and the dispersion relation u = Up implies u'/up = r>. 
However, the corresponding strength parameter a is found only after solving the full set 
of equations. As can be seen from this figure, there are no physical solutions for u' < Uq, 
two solutions for uq < u' < ui and only one physical solution for u' > Ui. (The full 
expressions for wq.i are cumbersome, but for /i > cr^/^ ^ 1 (corresponding to a mildly 
supermagnetosonic outflow), one find s ojn ^ ojp | l + (4/i^ — a^) a/ {S2fi^) + O (/i^^^'^)] and 



ui ^ Up ifi/sf^ [1 + (/i-'/')] , (cf. jiid 



2010l ).) Thus, the imposition of a finite energy 



flux per particle prohibits wave propagation at frequencies close to the cut-off u' = Up 
found from linear theory. At intermediate frequencies, waves of two different amplitudes, 
and different values of u^o are available to transport the required fluxes. One of these 
disappears (a 0) at finite u'/up, leaving only a single solution of the jump conditions at 
high frequency. 

In the case of a pulsar wind, the wave frequency is dictated by the rotation of the 
neutron star, but the plasma density decreases with increasing distance from the star. Thus, 
rather than solving the jump conditions as a function of frequency, more insight can be 
gained by solving as a function of radius. The connection between flux density and radius 
is found in terms of the luminosity per unit solid angle of the (radial) wind, dL/dQ^: 

{F') = \§- , (27) 
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which can be expressed in dimensionless form in terms of the quantity 



Ane^ dL 



(28) 



Physically, ol is the strength parameter of the circularly polarized vacuum wave that would 
be needed to carry the entire pulsar luminosity at the surface p = 1. For a spherically 
symmetric wind ol = 3.4 x lO^^Lgg^, where L38 = L/ (10^^ ergs^^). This quantity is related 



to the multiplicity parameter k, define d as the Goldreich 



cylinder divided by e (see, for example 



Lyubarsky fc Kirk 



ulian charge density 



2001 



at the light 



Kirk fc Mochol 



2011ah by 



K = ai,/ (4/i) 



(29) 



Writing the energy flux density as /x (J') leads to the expression 



{J') 



2nocT'i,u'^a1 



(30) 



Applying this to the case of circular polarization, where (J') = 2nQu'^Q and uj 



p' 



one finds u'^q = T^al/ {np^)- Then, the inequalities u'^q < l' < IJ^ and r> > 1 imply that 
circularly polarized waves propagate only when p > ai^/p. We therefore introduce the 
dimensionless scaled radius 



R 



PP/dL 

pn/A 



(31) 



which can be constructed after solving the jump conditions from the expression 

r>a; 



R 



2/i 



V 7o((-^')/^o7o) 
for both circularly and linearly polarized modes. 



(32) 



Figure [2] shows the properties of circularly polarized waves as functions of R for 
parameters appropriate to pulsar winds: a = 100 and p = 10100 (and, therefore, r< = 100). 
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Near the cut-off radius, which hes close to i? = 1, the mode properties are similar to those 
seen in Fig. [H except that the refractive index, now plotted function of radius, is 

no longer d egenerate. The detailed properties of these solutions have been discussed in 



KirkI ( I2OIOI ). However, note that these curves do not describe the ra dial evolutio n of a 



Asseo et al. 



f ll984h 



wave packet, as erroneously suggested in that paper and assumed by 
(who considered linear polarization), but simply specify the wave modes into which a 
wind of given /i and a can convert at a given radius. The difference arises because the 
radial momentum flux densit y, is not a conserved quantity in the radial evolution of the 
superluminal modes (see also 



Kirk fc Mochol 



2011bD. 



3.2. Linear polarization 

When the phase-averaged magnetic field vanishes, the properties of linearly polarized 
modes are similar to those of circularly polarized modes. This is illustrated in Fig. [31 where 
one sees that the refractive indices are practically identical. 

This figure also shows the nonlinearity parameter q. For linear polarization, q is defined 
in terms of quantities measured at phase zero. This point in the w ave is, however, specia l. 



197l|) 



It corresponds to the turning points of the "saw-tooth" waveform (iMax fc PerkinsI 
where the fluid velocity either vanishes or lies precisely in the propagation direction. These 
points do not play a large role in determining the average properties of the wave. As a 
result, q varies substantially, especially close to points where u^o changes sign (log(-R) ~ 0.4 
in the figure). One consequence is that an expansion of the fluxes using g as a small 
parameter is inadequate over most of the range relevant for pulsars, despite the fact that 
the waves are highly nonlinear (see the discussion in §1]). In the case of circular polarization, 
the quantities that enter into the definition of q (Eq. [7]) are phase-independent. The value 
of g is a useful characterization of the mode in this case, and Fig. [3] clearly illustrates that 
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one of the two solutions for given i? is a relatively weak wave, with q ~ few, whereas the 
other is stronger, with q ~ 10~^. 

The properties of linearly polarized waves which carry a non-zero phase-averaged 
magnetic field in the toroidal direction are shown in Fig. HI The refractive index is plotted 
for superluminal waves that correspond to /i = 10100, cr = 100, and four different values of 
the parameter 6, that describes the relative amount of phase-averaged magnetic flux carried 
in the subluminal wave (see Eq. |25|) . It can be seen that the cutoff moves to greater R as b 
rises. In the striped wind model, this corresponds to moving away from the equator, where 
6 = 0. 

Although the refractive indices for non-zero b have similar shape to those for 6 = 0, 
there is one potentially important difference: for b > 0.65 the lower branch cuts across the 
/3> = axis, and continues towards solutions with negative phase speed. Zero refractive 
index indicates that the H-frame coincides with the lab. frame. These solutions are standing 
waves in the pulsar wind, with constant magnetic field and wavelength ^ tlc- Negative 
refractive index indicates propagation towards the pulsar, in the sense that the velocity 
of the H-frame is inward propagating. Nevertheless these modes still carry the particle, 
energy, momentum and magnetic fluxes outwards towards the nebula. 

Conversion of a subluminal wave into a superluminal wave is accompanied by a 
substantial transfer of energy from the fields to the particles. This is quantified in Fig. [5l 
where the phase-averaged Lorentz factor (7') seen in the lab. frame is plotted as a function 
of R, for the same values of fi, a, and b as in Fig. |H For the smaller b values, 6 = and 
6 = 0.2, the upper branch of the solution has (7') ~ fi, which implies that almost all of 
the field energy is converted to particles in the equatorial regions of the wind. At higher 
latitudes, only a part of the Poynting flux can be converted, since the phase-averaged 
magnetic flux cannot be dissipated. However, even at 6 = 0.95, the upper branch of the 



.2 1—1 L, I 1 I 1 I I 1 

1 2 3 4 5 

log(R) 

Fig. 2. — The properties of circularly polarized waves for /i = 10100, a = 100, as a function 
of dimensionless radius in the pulsar wind. The inverse hyperbolic tangent of the refractive 
index ck'/u' = /3> is plotted, in order to bring out both the nonrelativistic regime, /3> ^ 1, 
near R = 1 and the relativistic regime, r> ^ 1, at 1. 




Fig. 3. — A comparison of the refractive indices for circularly and linearly polarized modes 
for fi = 10100, a = 100. The nonlinearity parameter q, defined in is also plotted for 
each mode. 
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dispersion relation carries particles with a mean Lorentz factor one order of magnitude 
larger than those of the incoming, striped wind solution. 



3.3. Latitude dependence 

According to the above figures, superluminal waves cannot propagate in the inner 
parts of a pulsar wind. The radius at which propagation becomes possible is close to 
R = 1 for vanishing phase- averaged magnetic field, and rises with this quantity. However, 
in the striped wind model, the luminosity of the wind, and, hence, the parameter ol is a 
function of latitude. In particular, the energy flux carried by the fields is proportional to 
sin^6', where 9 is the colatitude, (measured from the rotation axis). The flux carried by the 
particles, on the other hand, is expected to be much smaller, but its angular dependence 
is uncertain. Assuming, for simplicity, that the ratio a of these fluxes is independent of 
colatitude, one has 

ah = aL,eqSin6' (33) 

OL,eq = 4.2 X lO^^Lgg^, is the value of a-^ in the equatorial plane. 

The dependence of the average magnetic field carried by the striped wind on latitude 
depends on the inclination angle a between the magnetic and rotation axes. Two current 
sheets are present in the striped wind, located at phases 

(j) = ±0sheet + '2mr 

= ibarccos (cota cot6') + 2nn , (34) 



where n is an integer (IKirk et al. 



20021 ) . At ^ = atan(cota), the sheets merge, and vanish 



at higher latitudes (smaller colatitudes 6). At these sheets, the magnetic field reverses 
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Fig. 4. — The inverse hyperbolic tangent of the refractive index /3> for hnearly polarized 
waves that correspond to an incoming striped wind with fi = 10100, a = 100, and four 
different values of the phase- averaged magnetic field b, as defined in Eq. ( I25i) 
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Fig. 5. — The mean Lorentz factor of the outflow particles (7') measured in the lab. frame, 
for the linearly polarized superluminal wave solutions. Particles in a cold MHD outflow 
with the same energy and momentum flux have a Lorentz factor of 100 in this example 
(/i = 10100, a = 100). For the smaller b values, (7') ~ fi, implying a substantial fraction of 
the field energy has been transferred to the particles. 
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direction, so that 

(/"'/'sheet /-^ \ 

- / d(f)+ d(f)] /n 

= 1 — 2arccos (cota cot6') /tt (35) 

At those latitudes where there are no sheets, the striped wind model predicts an azimuthally 
symmetric magnetic field, i.e., no wave. 

From curves similar to those plotted in Fig. H] it is possible, for each value of 6, to 
identify the critical radius inside which the mode cannot propagate. This can then be 
interpolated to give a function Rcrip). Together with the definition R^^ = Pcr/^/c^L, this leads 
to 

^ = smeR,,[b{e,a)] . (36) 

^L,eq 

This surface is plotted in Fig. [6l for the case of the perpendicular rotator (magnetic 
inclination angle a = vr/2), and, in Fig. [TJ for various values of a. 

The location of the termination shock is determined by the balance between the 
momentum fiux carried by the wind and the external pressure. If we make the reasonable 



assumption that the latter is independent of 



shock is given for a ^ 1 by (ILyubarsky 



atitude, the approximate location pts{0) of the 



20021 . Eq. 9): 



dp, ^ ^ 



+pI = Ph^^'O (37) 

where po is a constant. The normalization of the radius of the termination shock surface 
depends on the external pressure and has to be calculated for each individual object. In 
Fig. E] three examples are shown, labeled by the ratio of the equatorial radius of the shock to 
that of the critical surface, Peq,ts/Peq,cr- It can be seen that if this ratio is significantly larger 
than unity, the cut-off surface falls almost entirely within the termination shock, except 
for a very small region close to the pole, where both surfaces have a cusp. A qualitatively 
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similar conclusion applies also for different magnetic inclination angles, as shown in Fig. [T) 
for Peq,ts/Peq,cr ^ & few, almost the whole critical surface falls within the termination shock 
radius, apart from a small region at latitudes close to the magnetic inclination angle a. The 
plots have been calculated using \i = 10100 and cr = 100, and show one quadrant of the 
poloidal plane. The 3D surface is generated by reflection in the p axis and rotation about 
the 2;-axis. 

When the cut-off surface falls inside the termination shock, superluminal waves can 
propagate in the wind. It is therefore possible, if the ratio Pcq,ts/Peq,cr > 1, that different 
modes coexist in the same wind but at different latitudes. The MHD description is more 
appropriate for the higher latitudes (close to a), and the superluminal wave solutions exist 
closer to the equatorial plane of the wind. If req,ts/''^cq,cr < 1, no superluminal modes can 
be supported by the wind at any latitude. In the extreme case where req,ts/''"eq,cr ^ 1, 
practically the whole wind can support these waves. 

The above calculations have been conducted for p ^ 10^. In this case the critical 
surface in the equatorial plane is located at p^r ~ 10'' for the Crab pulsar, whereas the 
termination shock is located at roughly pts ~ 10^. For the Vela pulsar, the critical radius 
lies at Per ~ 10^, which is also well wi thin the termination shock, located for this pulsar at 



Pts ~ 10^ (jKargaltsev fc Pavlov 



20081 ) . This value of p, corresponds to pair multiplicities 



of K = 2.5 X 10^ and k = 2.5 x 



conventional estimates fide Jager 



0^ for the Crab and Vela, respectively, much higher than 



20071 ). which would place the critical surface even closer 
to the pulsar. Thus, at least for young, energetic pulsars, Pcq,ts/Pcq,cr ^ 1 and superluminal 
waves can propagate in essentially the entire equatorial wind. 

However, the situation is different for pulsars which are members of binary systems, 
since the wind of the companion star can provide an obstacle capable of sustaining the 
termination shock relatively close to the pulsar. One example is the system PSR B1259-63, 
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containing a 48ms pulsar in an eccentric orbit about a B2e star. At periastron, the 
separation of the stars is p ~ 4.3 x 10^, while at apastron it becomes p ~ 6 x 10^. Very 
high-energy (TeV) gamma-r ays were predicted to emerge fro m the termination shock of 



this pulsar ( Kirk et al. 



Abdo et al 



2005 



2009 



19991 ) ■ and subsequently were detected (lAharonian et al. 
20111 ) at binary phases close to periastron. Assuming p ^ 10^ in this young 
pulsar, the critical radius at the equator is pcr ~ 2.4 x 10^, which is greater than periastron 
but smaller than apastron separation. Therefore, if the termination shock created by the 
interaction of the stellar wind with the pulsar wind lies close to the B2e star, superluminal 
modes may play a role in determining its properties near apastron, but not near periastron. 
It would, therefore, be interesting to search for an observable diagnostic of this transition. 



such as an orbital modulation of the radiation from the shock front. 



4. Discussion 



Nonlinear superluminal waves in electron-positron plasmas have been the subject of 
many papers over the last few decades. Nevertheless, their role in pulsar winds is still in 
need of clarification. One of the sources of this confusion is rooted in the use of the "stro ng 



wave limit" to sirnplify the treatment of linearly polarized waves ( iMax fc Perkins 



Kennel et al. 



1973 



1971 



Kennel fc Pellat 



19761 ). This corresponds to taking only the lowest order 



terms in the parameter q defined in ([7j). To this order, the particle flux in the H- frame 
vanishes - particles are said to be "locked" into the wave in this frame. Furthermore, the 
energy densiti es of particles and fields are in equipartit i on to lowest order in q. These 



properties led 



Kennel et al. 



(11973h and 



Kennel fc Pellat 



(119761 ) to assert that a strong wave 



would impose a rather severe upper limit on the particle flux from a pulsar. 

The formulation we present in ^and apply in [|3] avoids this limit. It demonstrates 
that highly nonlinear superluminal waves can carry a substantial particle flux, although the 



0.5 1 1.5 2 2.5 

P 



Fig. 6. — The surface pcr outside of which superluminal waves propagate is shown for a 
perpendicular rotator {a = tt/2). The pulsar is at the origin and z is the rotation axis. For 
comparison, the termination shock is also plotted for three different external pressures, such 
that its equatorial radius, normalized to that of pcr, is 1, 3 and 10. 




Fig. 7. — The critical surface for various inclination angles a between the magnetic and 
rotation axes, superimposed on the termination shock surface whose equatorial radius is 3 
times that of the critical surfaces. 
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minimum radius at which they can propagate moves outwards as the flux increases. These 
waves exhibit a wide range of q values (see Fig. [3]). Particles are not locked into the wave 
motion in the H-frame, and energy is not necessarily in equipartition. Physically, the reason 
for this apparent contradiction is that q does not adequately characterize the wave solution, 
since it is defined at a highly non-representative point in the phase of the wave; there exists 
a large family of nonlinear waves that have no counterpart with g = 0. 



The q uestion of t 



l e stab ility of superluminal waves has attracted considerable 



attention. lAsseo et al. 



( Il978l ) found that radiation reaction would damp these waves 



within a few wavelengths, for parameters appropriate to the Crab pulsar. However, this 
conclusion was based on the g = limit discussed abo ve, and may not a pply to large q 



modes. Again working analytically in the limit g = 0, 



Asseo et al. 



( 1l980l ) found rapidly 



growing Weibel-like instabilities, but noted that they are stabil ized when the 



magnetic field is sufficiently strong. For circular polarization. 



Lee fc Lerchel (11978 



P 



lase-averaged 



19791) 



also found rapidly growing instabilities when the plasma is locked into the wave in the 



H-frame. 
question. 



'article-in-cel 



Leboeuf et al. 



simulations have been used by several authors to investigate this 



(119821 ) injected a self-consistent wave and found it to be highly 



unstable when excited just above the cut-off 



at higher frequencies. 



Skjaeraasen et al. 



requency, but much less so when excited 



(l2005f ) injected a nonlinear wave into a uniform 



plasma, and again found a strong instability that could be partially stabilized in the 
presence of streaming motion. Although our results do not address this question directly, it 
can be seen from the dispersion relations presented in Figs. |2] and H] that the H-frame moves 
relativistically with respect to the pulsar (r> ^ 1) for R > 10. An instability that grows 
locally in the H-frame with growth-rate comparable to the wave frequency in that frame, 
propagates a distance of wavelengths in the pulsar reference frame in one e-folding time. 
Thus, although such instabilities may be present, their effect on the structure of the pulsar 
wind or the termination shock might not be dramatic. 
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In [|3]we show that superluminal waves are capable of carrying the particle, energy 
and magnetic fields thought to be transported by a pulsar wind when they reach the 
termination shock of an isolated pulsar, except for a range of latitudes around the (rotation) 
pole. In fact, at most latitudes, there is a large range of radii between pcr and pts in 
which propagation is possible. Our results do not address the question of where conversion 
from a striped wind to a superluminal wave might occur. However, since these waves are 
subsequently damped, it would seem likely that they could be sustained only within an 
extinction length of the termination shock itself, in which case their excitation and damping 
can be regarded as forming a part of this structure. 



LyubarskyI (j2003[ ) suggested the physics of the "termination shock" should differ 
from that of an MHD shock and permit reconnection of the alternating component of the 
incoming magnetic fi eld. A quantitative invest igation based on analytic considerations and 



ID PIC simulations (jPetri fc Lyubarsky 



20071 ) found that the effectiveness of reconnection 



depends on the ratio of the Larmor radii of electrons in the downstream flow to the 
wavelength of the stripes. Since this quantity increases monotonically with radius in a 
striped wind, these results imply that reconnection can be important only outside of a 
critical radius. Specifically, reconnection is predicted not to play a role for i? < 7r/12 (in 
our notation). At larger radii, the magnetic field is more effectively dissipated, until, at 
R > 7r-\/o'/3, full dissipation is achieved. However, this theory does not take account of the 
role of superluminal waves, although these waves are able to propagate in the region in 
which reconnection is predicted. 



Sironi fc Spitkovskyi (j201l[ ). on the other hand, performed 2D and 3D PIC simulations. 
They found that the shock triggered reconnection, leading to full dissipation of the magnetic 
field, over the entire parameter range investigated. In our notation, this range extends over 
1/20 < i? < 4. According to Fig. HI superluminal waves of very high phase speed (small /3>) 
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might be expected to play a role at the upper end of this range. However, the simulations 
do not appear to reveal their presence. The particle spectrum found in these simulations 
varies considerably, being closer to that observed in the Crab Nebula for lower values of R. 
However, as we show in ^ the termination shocks in the nebulae around isolated young 
pulsars are expected to lie much further from the pulsar than the range covered by these 
simulations, leaving us free to speculate that the physics controlling the structure and 
particle acceleration process might be substantially different in these objects. 

5. Conclusions 

We have reexamined the properties of transverse superluminal waves in an electron- 
positron plasma for parameters appropriate to pulsar winds. Our principal new results 
are: 

• The properties of circularly and linearly polarized modes with vanishing phase- 
averaged magnetic field are similar (Fig. |3]). 

• The widely used strong- wave approximation (g = 0), which underlies the assertion 
that these waves limit the particle flux, is inadequate. There is no wave-intrinsic limit 
on the particle flux. 

• The cut-off surface, inside of which these waves cannot propagate, has a strong 
latitude dependence. (Fig. [7]). It lies well inside the termination shock for young 
isolated pulsars, but must cross it at some phase of the binary orbit in the case of 
PSR B1259-63. 

• At high latitudes, where the phase-averaged magnetic flux carried by the wind is 
large, long wavelength, quasi- homogeneous modes of infinite phase speed (i.e., with 
[3 >= 0, Fig. Hj) can arise close to the cut-off surface. 
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On the basis of these results, and of recent discussions in the literature o f reconnection 



in the termination shock (jPetri fc Lyubarsky 



20071 : 



Sironi &: Spitkovsky 



201 ll ). we speculate 



that the physics of particle acceleration and the structure of this shock may be strongly 
influenced by the presence of superluminal waves. 



A. Two-fluid equations 



Following 



Clemmowl (Il977| ) and 



Kennel fc Pellat 



( 119761 ) ■ we present the solutions to 



the equations describing a cold, two-fluid plasma where transverse electromagnetic waves of 
superluminal phase speed can propagate. The equations are written in the H- frame and a 
plus (minus) index denotes positrons (electrons). It is convenient for the calculations to 
introduce complex quantities for the transverse four-velocity and the electric and magnetic 
fields: 



u±,± = Uy^± + iu^^± 
E = Ey + iE, 
B = By + iB, 



(Al) 
(A2) 
(A3) 



Using these, the continuity equation and the equations of motion for each species, expressed 
in the H-frame, are: 



d 
It 



{l±n±) = 



du\\ 



7 



7- 



dt 

du±^± 

dt 



Im(Mx,±5* 

mc 



mc 



-{'jE + iui±B) 



7^ = —Re{u^,^E*). 
dt mc 



(A4) 
(A5) 
(A6) 
(A7) 



where e is the magnitude of the electron charge. 
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In the H-frame, all space derivatives disappear. Conservation of magnetic flux, 
V ■ B = 0, is automatic, and the remaining Maxwell's equations (Ampere, Faraday and 
Coulomb law, respectively) become: 

1 dE 

-—— = — 47re(n+M_L^+ — n_M_L^_) (A8) 

M = " '^^) 

47rp = 47re(n+7+ - n_7_) = 0. (AlO) 



where p is the charge density. From (lAlOp . the charge density must vanish in the H-frame. 



In combination with the continuity equation, we get 

n+7+ = no7o = ^^-7- , 

where no, 70 are constants. The restriction to purely transverse electric fields (which is 
implied in the way Ampere's law is expressed above) means that the radial current must 
vanish: n^Ux+ = Ux-n_. In this case, the net force acting on the plasma in the transverse 
direction vanishes, so that Uy^z+ + 'Uy,z- is constant. In keeping with the approximation 
of toroidal fields (corresponding to the z-direction in the plane wave approximation) and 
radial propagation (x-direction), we set this constant to zero, so that Uy+ = —Uy^ and 
Uz+ = —Uz-- Because of these relations, we can restrict ourselves to solving the equations 
for the positively charged fluid, and will henceforth drop the plus subscript for simplicity. 
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A.l. Circular polcirization 

The simplest solution of these equations is the case of a circularly polarized wave. 
Then B — 0, and we have: 

u± = u±oe''^^* (All) 
E = Eoe^'^p* (A12) 
u\\ = uq (A13) 
n = no . (A14) 

The frequency coincides with the (proper) plasma frequency: 



Up - \ (A15 

V m 

and the magnitude of the perpendicular component of the four-momentum equals the 
strength or nonlinearity parameter a, defined using the electric field amplitude as 

u±o — a = — — — . (A16) 
mcujp 

The parallel component of the four- velocity is a constant. 

In the laboratory frame, which moves at speed —c/v^f, as seen from the H-frame, the 
frequency of the wave is 

oj' = r>a;po (A17) 

where r> is the Lorentz factor of the transformation between the H-frame and the lab. 
frame. 



A. 2. Linear polcirization 

Near the equatorial plane of a pulsar wind, the wave is linearly polarized and the 
relevant field components are the polar electric field and the toroidal magnetic field. In the 
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H-frame then, this requirement corresponds to Ez = and By = 0, i.e., E = Ey becomes 
real and B = iB^ purely imaginary (and constant). It follows that the z-component of the 
four velocity is constant lm{uj_). In a pulsar, a non-zero Uz corresponds to an azimuthal 
current loop, and, therefore, is proportional to the integrated magnetic flux crossing the 
enclosed surface. In the split-monopole geometry, this quantity falls off inversely with 
radius, and we set it to zero, so that u± is a real variable. 



Ampere's law becomes 

dE mccu u± 
d(f) e a7 

where 6 = ut and 



(A18) 



a = . (A19) 

Multiplying the equations of motion ( ]A5I) and ( ]A6ll by n and integrating, the four-velocity 
components can be expressed in terms of the variable y = E/Eq, the electric field normalized 
to the maximum value of its modulus (i.e., — 1 < y < 1) as follows: 

M|l = uo + aa^Bz{l - y) / Eo (A20) 
7 = 70 + ^ (1 - 2/^) (A21) 
u\ = '^'^ -u\-l . (A22) 

Having expressed all other unknowns in terms of the electric field, it remains to 
calculate the phase dependence of which is given by the equation 



« « ( — ) =^ (A23) 



or, explicitly: 
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( iKennel fc Pellatlll976[ ) with g, Q and h given by the expressions: 

470 

V 70 

-I- 

and 

N{y) = (1 - yf[{y + 1)^ - - g] + (1 - y)gQ 

is a quartic in y. Periodic waves exist if has four real roots, which are (in order): 

l/i = Ccos (e - 47r/3) - 1/3 
2/2 = C cos (e - 27r/3) - 1/3 

Z/3 = 1 

1/4 = Ccos(0 - 1/3 



(A25) 
(A26) 
(A27) 



(A28) 



(A29) 
(A30) 
(A31) 
(A32) 



where 



C = 2(4 + 12A2 + 3g)^^^ /3 
(16/27 + 2g/3 - 2AgMo/7o - 16AV3) 



^ = arccos 
provided ^ is reaL 



3^3 



^/2{^/?,u) 



3/2 



(A33) 
(A34) 



The nonhnear dispersion relation follows by demanding that the change in phase in a 
half-cycle of the field oscillation equal vr: 

dy 

ly, dy/d<p 



71 



(A35) 



The phase-averaged value of a quantity A that depends on y is given by 

1 f A{y) 



(A) 



7„i dy/ d(j) 



dy 



(A36) 
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If the phase-averaged magnetic field vanishes (A = 0), N{y) reduces to a quadratic 
in y^, with ?/2 = ~1 and = 1. In this case, the integrals in flA35p and (IA36P can be 
expressed in closed form in terms of elliptic integrals, at least for the functions of interest. 
For non-zero A, on the other hand, y oscillates between y2 > —1 and y^ = 1. The dispersion 
relation and the phase- averages must then be evaluated by numerical integration. 

The components of the stress energy tensor consists of a fiuid and a field contribution. 
In the H-frame these are: 

T°°, = 2mc'{n^'') = 2^070^0^(7) (A37) 
Tpart = 2mc^{n-fu^) = 2no-fomc^ (u^) (A38) 
T'l, = 2mc^nul) = 2no^omc'(^) (A39) 



7 



for the fiuids, and 



_ c{E)B _ c{y)\El 
^EM - - (A41) 

^EM = ^EM (A42) 

for the fields. For the circularly polarized modes, A = 0, {y) = 0, and (y^) = 1. In addition, 
Ux = and 7 = 70 are phase-independent. 
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B. Lorentz boosts 

The phase-averaged particle flux, energy, energy flux, and x-momentum flux densities 
and the electric field transform according to: 

J' - r> (2/3>noO'o + -/) (Bl) 

T'™ = r2^(r°° + /3>ro^ + /3^ri^) (B2) 

= r^^ [(i + /32)rii + ^>(r°° + r^^)] (B3) 

^,11 _ r^^ (T^i + 2/3>T°i + /32 rO°) (B4) 

{E') = r> (^>A£;o + (£;)) (bs) 
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